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Baryon number violating (BNV) processes are heavily constrained by experiments searching for 
nucleon decay and neutron-antineutron oscillations. If the baryon number violation occurs via the 
third generation quarks, however, we may be able to avoid the nucleon stability constraints, thus 
making such BNV interactions accessible at the LHC. In this paper we study a specific class of 
BNV extensions of the standard model (SM) involving diquark and leptoquark scalars. After an 
introduction to these models we study one promising extension in detail, being interested in particles 
with mass of 0{TeV). We calculate limits on the masses and couplings from neutron-antineutron 
oscillations and dineutron decay for couplings to first and third generation quarks. We explore the 
possible consequences of such a model on the matter-antimatter asymmetry. We shall see that for 
models which break the global baryon minus lepton number symmetry, (B — L), the most stringent 
constraints come from the need to preserve a matter-antimatter asymmetry. That is, the BNV 
interaction cannot be introduced if it would remove the matter-antimatter asymmetry independent 
of baryogenesis mechanism and temperature. Finally, we examine the phenomenology of such models 
at colliders such as the LHC. 

PACS numbers: 12.60.-i, ll.30.Fs, 14.80.Sv 



I. INTRODUCTION 

Experiments so far have not detected any violation of 
baryon number (_B), an accidental global symmetry of the 
standard model Lagrangian. Its conservation is dramati- 
cally demonstrated by the stringent limits on the proton 
lifetime, roughly r > lO'^^ years. The proton is not the 
lightest colour singlet, but it is the lightest colour singlet 
with baryon number, the conservation of which at low 
energy leads to its stability. 

Similarly stringent limits are set from searches of 
baryon number violating neutron decay modes, neutron- 
antineutron oscillations, and dinucleon decays. Baryon 
number violating processes must therefore be suppressed 
at low energy. The picture, however, may change at 
higher energy scales. There is no reason for B to be con- 
served; on the contrary, there are many strong theoretical 
reasons why this global symmetry should be broken. 

The SM already contains violation of baryon num- 
ber at high temperatures, through the non-perturbative 
sphaleron process, which breaks B and lepton number 
(L), but leaves {B — L) intact. There may exist copious 
other examples at high energy. Perhaps the most com- 
pelling theoretical reason for baryon number violation 
comes from the need to explain the matter-antimatter 
asymmetry. 

Observations of the cosmic microwave background 
and primordial deuterium abundance gives the ratio of 
baryon number density n;,, to photon number density 
of [II [5]: 

— = (6.1 ±0.3) X 10-1". (1) 
In a symmetric universe with _B = 0, the nucleons and 



antinucleons would continue annihilating until [3l |4] : 

^ = 10-19. (2) 

Assuming the universe starts off in a symmetric state, 
there must be a dynamical mechanism (baryogenesis), 
to create an asymmetry in baryon number. Baryoge- 
nesis presumably occurred at some high temperature, 
T > lOOGeV, and possibly much higher. 

Sakharov listed the three required conditions: (i) 
baryon number violation, (ii) C and CP violation, and 
(iii) departure from thermal equilibrium [5]. Though 
there is some room for the non-perturbative sphaleron 
process to be the baryogenesis mechanism (electroweak 
baryogenesis) in SUSY scenarios, the parameter space is 
highly constrained. Alternative scenarios for baryogen- 
esis have therefore also been proposed (for a review see 

ISIIII)- 

Furthermore, baryon number violation may arise quite 
generically in new high energy physics. In this case, there 
exists not only the BNV process which gives rise to baryo- 
genesis, but also other sources of BNV as well. Though 
not the mechanism for baryogenesis, these processes can 
affect the matter-antimatter asymmetry, even washing 
out an existing asymmetry [5]. We can use this to place 
additional constraints on new physics models or alterna- 
tively, if such additional BNV were detected, it may even 
be possible to rule out certain baryogenesis mechanisms 
as the source of the observed matter-antimatter asym- 
metry, as their contribution would simply be washed out 
by the additional BNV processes. 

Leptogenesis, for example, creates an asymmetry in 
the lepton sector at high temperature. This asymmetry 
clearly has a {B — L) component. The L asymmetry is 
then reprocessed into an asymmetry in the baryon sector 
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by the rapid {B + L) violating, but {B — L) conserving, 
sphaleron process. The introduction of additional BNV 
but {B — L) conserving interactions simply increases the 
rate of reprocessing. The addition of rapid (B — L) vio- 
lating interactions, however, would lead to a washout of 
the entire asymmetry. 

In this paper, we shall introduce baryon number vi- 
olation by first augmenting the SM with scalar di- 
quarks and leptoquarks. Introducing interactions be- 
tween these scalars then leads to baryon number vio- 
lation. Such interactions were originally motivated by 
charge quantization [9 , and have been studied as a pos- 
sible mechanism for barvogenesis|10j. 

Since the diquarks and leptoquarks all carry electro- 
magnetic and colour charge, the gauge interactions are 
too strong for the requisite departure from thermal equi- 
librium needed for baryogenesis to take place unless the 
new scalars have masses ^ O(TeV). In this paper we are 
interested in the phenomenology of these models at the 
LHC, that is with 0{TeV) masses, so we will not be con- 
sidering these interactions as a baryogenesis mechanism 
here. 

BNV interactions involving the first generation quarks 
are strongly constrained by nucleon stability. It is pos- 
sible to use the nucleon stability results to put stringent 
limits on dimension six BNV operators (the lowest di- 
mension effective BNV operators possible) involving the 
higher generations, meaning that BNV interactions are 
unlikely to be observed at the LHC [H] . If some unknown 
GIM-like cancellation mechanism were to suppress the 
nucleon instability, however, BNV interactions involving 
the top quark may become accessible at the LHC [T^. 

Instead of searching for a GIM-like cancellation mech- 
anism, by examining the flavour structure of new physics 
models, we concentrate on BNV that proceeds through 
multiple heavy scalars. Each additional scalar of mass M 
suppresses the amplitude of the nucleon decays by a fac- 
tor of Af^, so the nucleon stability experiments will put 
less stringent conditions on M. We shall be looking at 
specific models, extending the SM with renormalizable 
terms. (The effective operators would then correspond 
to dimension nine or twelve, instead of the dimension six 
operators considered in Refs. fTTl [T^ .) 

In Section|ll]we review the possible scalar leptoquarks 
and diquarks and their BNV interactions. We examine 
the constraints from nucleon stability on one of the in- 
teractions, considering couplings to first and third gener- 
ation quarks. In Section HI we analyse the effect such an 



interaction can have on baryogenesis, showing stringent 
limits apply for models which break the global {B — L) 
symmetry. Finally, in Section |lV| we examine the phe- 
nomenology of this interaction at the LHC. Detection of 
such BNV interactions may then allow us to rule out cer- 
tain baryogenesis mechanisms, even if their mass scale 
lies outside the reach of present day colliders. 



II. BARYON NUMBER VIOLATING SCALARS 

A. A catalogue of models 

The models studied in this paper are those introduced 
in Ref. and were motivated by charge quantization. 
Here we use the same notation for the list of all possi- 
ble scalar leptoquarks and diquarks, and their transfor- 
mation under the SM gauge group SU{'i)c x SU{2)l x 
U{1)y- The convention here is that the scalar fields have 
the same quantum numbers as the quark and lepton bi- 
linears: 



CTi.i ~ QhifLY ^ uiiicRY ^ dRivity 

^1.2 ^QlUlY 

0-3.1 ~ QjAQlY ^ u]^idRY 

^3.2 --QiAQlY 

fa. 3 ~ QjAQlY ~ URidnY 

03A -^QiXQlY 

CTe.i ^ U]^{urY 

(76.2 ^ ur{urY 

0-7.1 ^ diiidRY 

07.2 ^ dnidnY 

(Tg ^ dnieiiY 



(3,1,2/3) 
(3,3,2/3) 
(3,2,-7/3) 
-2/3) 
-2/3) 
-2/3) 
-2/3) 
-4/3) 
-1/3) 
-8/3) 
-8/3) 



(3,1, 
(3,3, 
(6,1, 
(6,3, 
(3,1, 
(3,2, 
(3,1, 
(6,1, 



(3,1,4/3) 
(6,1,4/3) 
(3,1,8/3), 



(3) 

where ^E"^ = C^*"^, and C is the charge conjugation ma- 
trix. The standard model fermions and right-handed neu- 
trinos transform in the usual way: 

/i^ (1,2,-1), e,j~ (1,1,-2), z/fl^ (1,1,0), 
- (3,2,1/3), u;^- (3,1,4/3), - (3, 1, -2/3). 

(4) 

Note the leptoquarks in Eq.([3]) carry B = —1/3, and the 
diquarks B — —2/S. As we are breaking this symme- 
try, the following particles then carry identical quantum 
numbers: 



0-1.1 

0-1.2 



'3.1 
''3.2 



0-4 = <7j^ 

O6.1 = o-g 



-(3,1,2/3) 
-(3,3,2/3) 
(3,1,-4/3) 
(3,1,-8/3). 



(5) 



Such particles with both leptoquark and diquark cou- 
plings leads to baryon number violation. For example, 
taking the a = cti.i = (J^ i ^ (3,1,2/3) model the La- 
grangian is extended to include the terms: 

C D Ai.i(e_R)=a-MH 4- X3.i{uR)(7d'p. + H.c, (6) 

where A 1.1 and A3.1 are dimensionless coupling con- 
stants, and generational and colour indices have been 
suppressed. If we take both Yukawa couplings to the 
SM fermions to be A = Ai.i ~ A3.1, we obtain the fol- 
lowing estimate of the proton decay rate on dimensional 
grounds (see fig IT]) [9]: 



r - O 



A^MpS 
Af4 



(7) 



3 



where Mp is the proton mass. Given the bound on the 
partial Ufetime for p — ^ 7r°e+, r > 8.2 x lO'^'^ years [l3] . 



this translates to a limit on the mass of: 
M„ > A X lO^^GeV. 



(8) 



will not investigate the possibility of such a mechanism 
here. Instead we restrict our attention to the scalars 
which are not subject to this large mass bound, namely: 

CTS.S, Cr3.4, fS, a'6.2, Cr7.2. 

Baryon number is then broken by the introduction of 
two scalars and an interaction between them. The pos- 
sible scalar-scalar interactions, without the particles of 
Eq.(|5]), are listed below. The AS — 1 list is: 





a 



(12) 



where 0" is the physical SM Higgs scalar, p represents 
a new Higgs hke scalar p ~ (8,2,1), and crfjCrf are the 
components of the 0-5 SU (2) ^ doublet. The AS — 2 list 
is: 



FIG. 1. Proton decay, p — > vr e , in a one exotic scalar ex- 
tension with coupling to the first generation. 



''■3.30'3.30'7.2 
O'3.40'3.40'7.2 
O'6.20'7.20'7.2. 



(13) 




FIG. 2. Proton decay, p ^ K+ V, in a one exotic scalar ex- 
tension with coupling to the third generation. 

If we instead take the dominant coupling to be to the 
third generation, proton decay now proceeds through the 
diagram in Figj2j The decay rate is estimated as (see 
appendix for details): 

V4„8 



Mi 



\Vub\\Vtd\\VtsmM, 



In 



My, 



In 



Mt_ 

Mr 



(9) 



where is the weak coupling constant. Given the 
limit on the partial lifetime for p — K^v, t > 2.3 x 
lO'^'^yearsjiMj, this translates to a limit on the mass of: 



M„ > A X 10"GeV. 



(10) 



This corresponds to "opening up" the effective operator 
O^ull - {Jb^'tR){J^'lR) (11) 



of Ref.[TT] and is comparable to their limit on its coef- 
ficient C[% < 10""TcV"^. Such BNV physics is then 
only accessible at the LHC if we assume a GIM-like can- 
cellation mechanism to avoid these stringent bounds. We 



Note that the models of Eq.(12) conserve (B — L) while 
the models of Eq.( 13 ) break [B—L). This is an important 
distinction when it comes to consideration of the interac- 
tion's effects on baryogenesis. Breaking baryon number 
using these scalar-scalar interactions leads to less strin- 
gent bounds from nucleon stability. 

Considering this shorter list, we will study the 
"■3. 3173. 3177. 2 interaction in greater detail, due to its in- 
teresting phenomenology. First of all, it involves di- 
quarks, which can be produced more favourably than 
leptoquarks in hadron colliders. Secondly, (T3.3 couples 
to an up and down type quark. This means the possi- 
bility of top quarks in the final state, which is a clearer 
signal than generic jets. Finally, for simplicity, we ex- 
amine the SU(2)]^ singlet, tT3.3, rather than the SU{2)]^ 
triplet, (73.4. 

The a particles in Eq.(|3]), are either scalar diquarks or 
scalar leptoquarks. The usual constraints, from low en- 
ergy experiments, on the couplings of diquarks and lep- 
toquarks to standard model fermions then apply. For 
example, from consideration of atomic parity violation 
of caesium, leptoquarks coupled to the right handed first 
generation fermions have a bound M > A x 2TeV |15j . 
For further constraints on the leptoquarks (in particu- 
lar on products of different couplings) see Refs.[15 18j. 
Products of diquark couplings are also constrained from 
limits on flavour changing neutral currents (FCNC) and 
meson-antimeson oscillations |19) . These limits, however, 
do not apply in the limit where all but one of the cou- 
plings vanish. 

There are also constraints on the masses of leptoquarks 
and diquarks from collider experiments. From searches 
of pair produced leptoquarks the best limits for the three 
generations are 



(14) 



Mlq > 384GeV 1st generation {jjee) 
Mlq > 420GeV 2nd generation (jj>/x) 
Mlq > 148GeV 3rd generation (56z/i/), 
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where we have indicated the final state searched for with 
j =jet. 

Diquarks can be produced singly at hadron colliders 
and then decay to two jets. It is then possible to exclude 
a given diquark given its coupling and mass from searches 
of dijet resonances [24l - f26] . Unlike the searches for pair 
produced leptoquarks, however, this diquark search is 
model dependent as the production occurs through the 
Yukawa coupling to the SM quarks. 

As noted before, such baryon number violating scalar 
interactions can also be the mechanism for baryogenesis 
[TU] , The (73. aCTs. 3(77. 2 model, studied in the rest of the 
paper, was proposed as a baryogenesis mechanism by the 
authors of Ref . [27l [28] . The addition of a singlet scalar 
field, which decays out of equilibrium, allows the third 
Sakharov condition to be met, even with ©(TeV) masses 
for the coloured scalars. 

Alternatively, colour triplet scalars with both diquark 
couplings to the SM and leptoquark coupling to a SM 
quark and a right handed Majorana neutrino would also 
allow the third Sakharov condition to be met [29] . These 
papers also discuss the neutron-antineutron oscillations 
induced by such fields [50] . 

Here we do not consider these particles as a baryo- 
genesis mechanism, instead we obtain constraints from 
washout of baryogenesis. 

B. The Particular Model: (73.3(73.3(77.2 

In the remainder of this paper we shall focus on the 
(73. 3(73. 3(77. 2 model. This model has the advantage of not 
containing any of the conjugate pairs of Eq.([5|: recall 
that if the particle has both diquark and leptoquark cou- 
plings, the stringent limit on the mass of Eq.(|lO|) ap- 
plies. In addition the interaction (73.3(73.3(77.2 can produce 
top quarks in the final state which would give an inter- 
esting signal at the LHC. Furthermore, the production 
of single diquarks through the Yukawa coupling is more 
favourable than leptoquarks at hadron colliders. The SM 
Lagrangian is extended to include the terms: 

C D A3(dfl)=WflCr3. 3-1- A7((iij)'=(ij^(77. 2 + ^0-3. 30-3. 3f77.2 + ^^-<3- 

(15) 

where A3, A7 are dimensionless coupling constants and fi 
has dimension of mass. Generational and colour indices 
have been suppressed. The exotic fields transform in the 
following way under the SM gauge interactions: 

(73.3^(6,1,-2/3), (77.2^(6,1,4/3). (16) 

We will primarily be interested in the Ma-y > ^M^^ hier- 
archy, especially for phenomenology at the LHC, which 
we study in Sec|IV| 

The couplings A3 and A7 are constrained from low en- 
ergy experiments measuring FCNC and neutral meson 
oscillations. For example if the coupling to the third 
generation quarks A^^ = 0.1, the coupling to the second 
generation quarks is constrained to be Af^ < 10""^, for 



Mcr^ = ITeV or the — oscillation frequency would 
exceed the experimental value [19l |3T1 [32] . 



d . . d 
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FIG. 3. Tree level n — n oscillation for the (770^0^ model with 
the scalars coupled to the first generation quarks. 




FIG. 4. Two loop nn — >■ tt tt decay for the (77(73(73 model 
with the scalars coupled to the third generation quarks. 

For couplings to the first generation quarks, there is 
also a stringent constraint from neutron-antineutron (n— 
n) oscillations (see Figj3|. Using dimensional analysis, 
the rate of oscillation is estimated as: 

where M„ is the mass of the neutron. The lower limit 
from experiment for the oscillation time is currently 
around r > lO^s [Ml [31 ■ Taking A = A3 « A7, 
M„ = ^ Mc-j, and /i ^ \M„, we obtain a bound: 

M<T > A"*/^ X lO'^GeV. (18) 

A similar constraint from dinucleon decay has previously 
been calculated [S]. 

We have set ^ ~ XM„ as this ensures the two decay 
modes of a^ have similar branching ratio. In the limit 
/i — )■ 0, the BNV aspect of this model is switched off and 
no constraints from nucleon decay or washout exist. A 
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FIG. 6. Two loop n — n oscillation for the aja^az model with 
CT7 coupled to the first and 0-3 to the third generation quarks. 



III. WASHOUT OF BARYOGENESIS 



A. High Temperature Baryogenesis 



FIG. 5. Two loop nn — >■ tt^tt^ decay for the a-juzo-i model 
with the scalars coupled to the third generation quarks. 



small /i, however, makes the branching ratio for cfy — > 
C3 + (73 insignificant. Detecting any BNV interaction at 
colliders then becomes even more unlikely. 

The constraint on however, changes when the 

dominant coupling is taken to be to the third genera- 
tion SM fermions. The most str inge nt constraint then 
comes from nn — > vt^tt" (see Figs 4l5 1, which has a par- 
tial lifetime r > 3.4 x lO'^" vears |35]. Evaluating the 
integrals around the loops, and then using dimensional 
analysis and taking ^ ^ AA/^r, the double neutron decay 
rate is estimated to be: 



O 



This translates into a limit on the mass of: 
Ma > A'*/^ X lOGeV. 



(19) 



(20) 



The simplest n — n diagram now involves four loops, so 
we expect it to be even more suppressed. If the particles 
couple predominantly to the third generation, we see that 
nucleon stability constraints do not preclude BNV from 
occurring at LHC energies. We will soon see, however, 
that far more stringent constraints arise from washout of 
baryogenesis. 

Finally, let us mention an intermediate case where the 
(T7 couples to the first generation and the (T3 to the third. 
The limit on the mass again comes from n — n oscillations 
(see Figj6]). The oscillation rate is estimated as: 



r(n-n) - O /xAyA^ 



Mt 



M2 M4 
(T7 0-3 



This translates as a limit on the mass of: 
Ma > A*/^ X 6TeV. 



If baryogenesis were to occur before the temperature, 
T, of the universe reached the mass scale of the a par- 
ticles, as would be the case in some high temperature 
leptogenesis scenario, the creation and decay of the a par- 
ticles in a A(i3 — L) 7^ sequence can erase the matter- 
antimatter asymmetry. For example, in the Ma-, > 2Ma^ 
case such a sequence could proceed on-shell in the follow- 
ing way: 



d + — > (77 

0> ^ (T3 -I- CT3 

(T3 — > u -I- d. 



(23) 
(24) 
(25) 



If the rate of one of these steps, F, is less than the expan- 
sion rate of the universe, H, the A{B—L) ^ process will 
not be occurring rapidly enough for washout [H [37] . 
If we again take the A = A3 ~ A7, /i ~ XMa case we see 
to avoid washout we require: 



A^M 



Ma 
T 



9 



1/2^2 



M^ 



PI 



(26) 



where g specifies the degrees of freedom, Mpi is the 
Planck mass, and we have included a Lorentz factor for 
the typical rate. Remembering that this is required to 
hold for all T > Ma (below which the initial inverse de- 
cay will be Boltzmann suppressed) , this translates into a 
washout avoidance condition on A: 



A < 



'/'Ma 



Mpi 



1/2 



(27) 



This is a more stringent constraint than Eq.(18l for 



©(TeV) masses. So washout of high temperature lepto- 
genesis could take place with the addition of such fields, 
even with couplings only to the first generation (which 
are the most constrained from nucleon stability). 

The addition of such fields with couplings greater than 
the above constraint, would mean any existing asymme- 
try in the quark sector above the mass scale of the a 
particles will be removed by their BNV interactions. 
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B. Low Temperature Baryogenesis 

Now we examine the scenario where baryogenesis oc- 
curs at a temperature below the mass scale of the a par- 
ticles, which may be the case e.g. in electroweak baryo- 
genesis. The inverse decay dd — > <j-j is now Boltzmann 
suppressed: 



Ma 

T 



3/2 



exp 



T 



(28) 



To avoid washout we require this rate to be less than H 
at the temperature of baryogenesis, T;,, giving the bound: 



A < 



exp 



n 



9 



1/2 



MpiM, 



5/2 



(29) 



Setting Mo- 
bound of: 



= ITeV and T,, = lOOGeV, we obtain a 



A < 10" 



(30) 



Due to there being only one order of magnitude differ- 
ence between M„ and T;,, the inverse decay is not yet 
sufficiently suppressed for the on-shell sequence to be in- 
significant, and thus the stringent bound on A, Eq.(29), 
applies. 




(77 




^3, 




FIG. 7. BNV scattering for the (770303 model 

If M„ were increased substantially (to « 4.5TeV), the 
inverse decay does become sufficiently suppressed. The 
off-shell BNV scattering process depicted in FigjT] then 
becomes dominant, with rate: 



r 



Ml 



(31) 



When compared to H, and again taking fi ~ A^M^, this 
translates into the washout avoidance condition: 



A < 



1/2^1 



T^Mpi 



1/8 



(32) 



which is more stringent for = lOOGeV than the con- 
straint from Eq.([29| only for M^ > 4.5TeV. For these 
larger masses, we can have A ^ 0.1 — 1 while still sat- 
isfying constraints from washout avoidance and nucleon 
stability. 

For lower masses M^ < 4.5TeV, the couplings of this 
{B — L) violating interaction are tightly constrained by 
the washout avoidance condition, no matter to which 
generation quarks the scalars couple. A stronger cou- 
pling would mean baryogenesis would be washed out if 
it occurred at a temperature > lOOGeV. As we are 
almost certain it occurred at a temperature of at least 
IGOGeV, we take the constraint of Eq.(30l as our start- 



ing point for our examination of this model at the LHC 
for masses below 4.5TeV. We will return to the high mass 
case. Mo- > 4.5TeV, at the end of the next section. 

Finally we point out that if A lies between the high 
and low temperature washout constraints, detection of 
such interactions could be used to rule out certain high 
temperature baryogenesis mechanisms, without having 
to probe their possibly high energy scales ~ lO^^GeV 
directly. 



IV. COLLIDER SEARCHES 

The production of the exotic scalars can take place ei- 
ther through the Yukawa coupling with the quark fields 
or by gluon-gluon fusion. Let us first examine both pos- 
sibilities, in the case where the heaviest diquark mass is 
Mar ^ 4.5TeV, in light of the stringent washout avoid- 
ance constraint in Eq. ( 30 1 . 
In the limit of /i 



0, we effectively turn off the 
BNV interaction and return to the standard diquark phe- 
nomenology. The Yukawa coupling to the first generation 
may then be large enough for the — cr — > jj signal to 
be significant. Such diquark production has been stud- 
ied in [311 I3H1 ES]- Experimental limits on such dijet 
resonances have already been set P^26j . 

If /i ^ AA/o-, the washout avoidance condition implies 
a lifetime cr > lO^^m, so the observation of displaced 
vertices may occur. The scalars would be pair produced 
through gluon-gluon fusion (see Figjs]), which depends 
only on the colour charge. For a unity branching ratio 
into two jets, the LHC at = 14TeV with C = lOOfb"^ 
of data, can discover such pair produced colour sextet 
scalars with masses below 1050GeV [301 HI]- If we as- 
sume the hierarchy M^^ > 2Ma-^ , we must only concern 
ourselves with the branching ratio for cfy —>■ (73 + CT3 and 
to which generation quarks the a particles decay. If the 
coupling was to the first two generations, it would not be 
possible to extract the existence of a BNV process, due 
to the generic nature of the jets. 

Reference [32] provides a plot of the pair production 
cross section of colour sextet scalars as a function of the 
particle mass. Pair production of the scalars declines 
from w 6pb for M = 500GeV to w 60fb for M = ITeV 
at = 14TeV. Taking Mc^ > ^M^^ , branching ratio r 
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for (J^ — > (T3(T3, and coupling of the scalars to the third 
generation quarks we have the final states of pair pro- 
duced 0-7177 in Table. I (also see Fig|8]for an example of 
the decay chain). 



Final state 


Fraction of total 


hbhh 




bb tb tb 


r(l — r) 


tb tb bb 


r(l — r) 


tb tb tb tb 





TABLE I. Final states from pair production of 0-7(77 given 
branching ratio r for cr? — >■ 0-30-3, and coupling of the scalars 
to the third generation quarks. 



from -washout: Eq.(32 ). It saturates the nucleon stability 
limit, Eq.(22), to the level of precision -with -which that 
limit "was calculated. 

Using the Breit-Wigner approximation for scattering 
through resonances and the CTIO parton distribution 
function [43], -we obtain the cross sections for tbtb and 
jj final states in Table. II (also see Fig|9|. 



M^7 (TeV) 


4.75 


5.0 


5.25 


o-(^ = 14TeV) (pb) 


0.19 


0.17 


0.16 


aiy^ = 7TeV) (pb) 


0.04 


0.04 


0.03 



TABLE II. Final states a{ibtb) = a-(jj) from 07 production 
for pp collisions at ^ = 14TeV and ^/s — 7TeV. The Yukawa 
coupling to the quarks is A7 = 1, and branching ratio 07 — > 
030-3 equal to 0.5. 






FIG. 8. Pair production of 0707 through gluon-gluon fusion. 
The 07 decays into 0303 resulting in a tb tb bb final state. The 
bb tb tb final state from the 0-707 pair production would occur 
in the same amount. 

One could then utilise the semi leptonic decay to dis- 
tinguish top from anti-top. This -would reduce the signal 
by a factor of (2/9) for every top quark required to decay 
to a lepton. But we are still left -with many b jets, the 
charge of "which is even more difficult to extract. Further- 
more due to the initial B = ajdj state, there are the 
same number of quark and antiquark jets over a many 
event average. Such a BNV signal is therefore unlikely 
to be detected conclusively at the LHC. 

Now to the case A — )■ and /i ^ XM^. Here the 
phenomenology does not change much, except for large 
enough /i, the decay width 07 0-3 -I- 0-3 becomes large 
and u^ may then not decay at a displaced vertex. Given 
the mass hierarchy we have chosen, M^^ > 2Ma^ , 0-3 will 
still decay preferably to two jets rather than through the 
virtual cry propagator. 

Finally we examine the case where the heaviest diquark 
mass is large: > 4.5TeV. The couplings are then far 
less constrained from washout and there may be signifi- 
cant production — > 07, as well as a significant branch- 
ing fraction dj (73 + (T3. For purposes of illustration 
we choose values A7 = 1 and r{(7j 0-30-3) = 0.5. Fur- 
thermore we assume 0-3 decays to the third generation 
quarks. The choice A7 = 1 is consistent with the limit 



Given that the limit on dijet resonances at M — 
4.0TeV at ^ = 7TeV is 0.005pb [25l ESI, such a high 
coupling A may very quickly be ruled out once more data 
is released (the experimental papers do not supply any 
limits for dijet resonances above 4.1TeV, we can imag- 
ine such a large dijet excess at 4.75-5.25TeV presumably 
would have already been pointed out if it existed) . This 
scenario should therefore be taken as an illustration of 
the BNV signals accessible once the centre-of-mass en- 
ergy is increased to y/s — 14TeV. Such a clear signal, 
however, relies on careful selection of the parameters and 
should therefore be taken as a best case scenario. 





FIG. 9. BNV at the LHC dd ^ 07 ^ 03 + 03. 



V. DISCUSSION 

Let us summarize the overall phenomenological picture 
for such BNV models. First the low mass regime for the 
heaviest scalar: A/g.^ < 4.5TeV. For the diquark inter- 



actions of Eq.(13), [B — L) is broken so for such a low 
mass, one of the couplings must be small for washout of 
baryogenesis to be avoided. This forces us to rely on pair 
production for any sizable BNV process to be occurring 
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at the LHC, making detection of the BNV interaction 
very difficult. Even without a clear BNV signal, multiple 
scalars may be discovered through the four jet pair pro- 
duction signal, or in the case — >■ 0, the dijet signal if A 
is large enough. 

Once the properties of the scalars begin to emerge, 
one would infer the possibility of a gauge-invariant cubic 
interaction between the recently discovered scalars. Thus 
the possibility of BNV in this sector of particles would 
be hinted at. The precision study of such an interaction 
would be left to further experimental work. 

If the heaviest scalar mass is large, M^^ > 4.5TeV, the 
washout condition is much weaker and a clearer BNV 
signal may occur albeit for a very specific choice of pa- 
rameters. If such a choice is not realised we are again 
left with signals of multiple scalars. For masses greater 
than around ITeV, only the dijet signal governed by the 
Yukawa coupling to first generation quarks is accessible 
at the LHC, pair production being too small for such 
masses. 

If we instead examine the {B — L) conserving inter- 



actions of Eq.(12) we are left with leptoquark instead 
of diquark scalars. The washout constraint then does 
not apply, but being leptoquarks, these again would not 
be produced singly in large numbers at hadron collid- 
ers (pair production through gluon-gluon fusion is still 
possible). That is, even in an optimistic scenario with 
EM strength Yukawa coupling, the production mecha- 
nisms g + d d* — > (T5b + e~ and pp^j + d + X^ 
e~^e~ +d+X — > +a5b+X , only result in a cross section 
of 3.4fb[li] and 10fb[Sl SS] respectively for M„ = ITeV 
at -^/s = 14TeV at the LHC. But due to the nucleon sta- 
bility constraints this coupling is expected to be small. 
So we expect BNV leptoquarks to only be pair produced 
in large numbers in a hadron collider (not produced on 
their own). This again results in a i? = initial state, 
making detection of the BNV interaction difficult. 



VI. CONCLUSION 

We have examined an extension of the standard model 
involving diquarks with renormalizable terms. We have 
shown that the nucleon stability constraints are greatly 
weakened if the exotic scalars are coupled to the third 
generation quarks. The diquark-diquark interaction vi- 
olates {B ~ L), so for a mass of the heaviest diquark; 

^ 4.5TeV, the most stringent constraint now comes 
from washout of baryogenesis. For such a choice of mass 
we are left with very stringent constraints on at least one 
of the couplings no matter to which generation quarks 
the scalars couple. Although the new scalars can be pro- 
duced at the LHC, these stringent coupling constraints 
greatly inhibit any clear identification of such BNV in- 
teractions at the LHC. If the heaviest diquark has a mass 
higher than 4.5TeV, the most stringent constraint again 
comes from nucleon stability. If the scalar coupling to 
the SM quarks then involves the higher generations, the 



couplings may be quite large, allowing for certain BNV 
processes to occur rapidly at LHC energies. In either 
case, even if no BNV interaction is detected conclusively, 
the exotic particles would show up in other ways in the 
experiments. The discovery of multiple scalars would be 
interesting in itself, and once an overall picture emerged, 
the possibility of BNV in this sector of particles could be 
postulated. 
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APPENDIX 

To estimate the rate of the process in Fig[2] wc ffi-st 
evaluate the contribution from one of the fermion lines 
forming a loop with the W boson. Using relations for 
the spinors such as m = Cv"^ (see appendix G.4 in 47'), 
and ignoring the external momenta we find the Feynman 
amplitude for such a loop: 



xAciR 



1: 

fc2 



nit 



g^,VtdY'Lu, (33) 



where L(R), is the left (right) chiral projection operator. 
This can be simplified to: 



(4A3VidKb5^mbmt) X (vLu) 
d^k 1 1 



1 



(27r)^ fc2 — — fc^ — 

Taking the leading term of the integral: 



( -) 



(34) 



(35) 



we see this fermion line contributes a factor to the overall 
amplitude of: 



1 rrib 

2tT^ TTlt 



In 



mt 



vLu. 



(36) 



Then ignoring the external momenta, we find the depen- 
dence on the propagator masses for the overall diagram. 
The matrix element is simply estimated to be the mass 
of the neutron, with exponent set to restore the correct 
dimension of mass to the decay rate. This yields the rate 
of Eq.(|9f. 
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To estimate the rate of the process in FigjS] we note 
the integral around the loops has the form: 



d^a / 1 



1 1 



d^fc 1 



1 1 



(27r)4 fc2 _ fc2 _ ^2 fc2 



2 — m?. }) 



(37) 



Focusing first on the integral over fc, and introducing the 
Feynman parameters Xi, we rewrite this integral in the 
following form: 



4:1 / / da;ida;2da;3dx4 

Jo Jo Jo Jo 



X / d^k 



|(fc + axi)^ + a^xi{l — xi) 
{1 — Xi — X2 — X:} — X4^)mf 

} — 5 

(38) 



(We have reintroduced the ze, e > term in the denom- 
inator of the propagators.) After changing variable to 
I = k + axi this becomes a standard integral which eval- 
uates to: 



-2i 



1 nl — Xi nl — X2 — Xi nl — X3 — X2 — xl 



(47r)2 Jq jq 

X |a^xi(l — xi) — (I ~ xi — X2 — X3 — Xi)m'f 



da:ida;2da;3da::4 

2 
t 

-3 



ximl^ - X2'm'^ - x^ml^ - Xiinl + iej 



(39) 



Now writing — UQ^af, the denominator is zero for: 

oo = ± I — T-^ rN?(l - xi - X2 - x^ ~ X4) 

lxi(l - Xi) 

1 1/2 

-I- ximl^ + X2m^ -f x^ml^ + Xiinl] +af> ^ ie' . 

(40) 



This means any singularities occur in the second or fourth 
quadrant on the ao plane, allowing us to Wick rotate so 
that ueq = iao and o^q = —a^. 
Equation (37 1 now becomes: 



-2 r d^as / l_ 

(47r)2 J (27r)4 \a% +ml aj^ + a% + mlaE + ml^ 

fl rl — Xi /•i — X2—Xi /•1 — X3 — X2—X1 

X I I I I dxida;2da;3da;4 



1 1 

,2 „2 ^ ^2 „2 _L ^2 



1 



"'0 











x{a|;a;i(l — xi) + {1 — xi — X2 — x^ — X4)mf 

+ ximl^ + X2m^ + xsml^ + X4,ml}^^^ . 

(41) 

Now this is in the form J daf{a)g{a), where g{a) rep- 
resents the integral over the Feynman parameters, with 
f{a),g{a) positive definite. The integral over the Feyn- 
man parameters, g{a), is maximal for a = 0. If 
g(a) < A for all a, and f{a),g{a) are positive definite 
J daf{a)g{a) < Jdaf{a)A. Applying this to Eq.( 41), we 
find the leading term of this integral to be: 



(42) 



where the second bracket is the leading term oi g{a = 0), 
and the first the remaining integral J da/ (a), and we 
have ignored factors of i, n, and logarithms of 0(1). Us- 
ing a similar technique of dimensional analysis as outlined 
before we again obtain the rate of Eq.( 19 ). 
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